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Abstract 



| A model in which the soft collisions of the nucleon are described in terms 

of interactions of its two constituents (a quark and a diquark) is proposed. 
When adjusted to describe precisely the elastic proton-proton scattering data 
and supplemented with the idea of wounded constituents, the model accounts 
rather well for the centrality dependence of particle production in the central 
rapidity region at RHIC energies. 

1 Introduction 

> ■ 

The data on production of particles in relativistic heavy ion collisions collected dur- 
ing the operation of the RHIC accelerator [1-3] indicate that the wounded nucleon 
model [4], does not describe correctly the observed particle multiplicities. Con- 
trary to the model predictions, the particle density per one wounded nucleon in the 
central rapidity region (i) increases slowly with centrality of the collision and (ii) 
substantially exceeds that observed in nucleon-nucleon collisions. 

This is not too surprising because the model was always considered only as a 
first, rough approximation. Indeed, soon after the original model was formulated, 
a possible improvement was proposed in the form of the wounded quark model 
[5]. Recently, it was suggested [6] that the wounded quark model can account for 
the main features of the data in the central rapidity region. To obtain this result, 
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however, it was necessary to assume a rather large number of quark-quark collisions 
in the nucleon-nucleon scattering which is difficult to justify. If more realistic value of 
this number is used, the model predicts larger multiplicities than actually observed 
[7]. Nevertheless, these results show that the idea of a wounded constituent model 
may be not far from reality. 

In the present paper we propose another generalization of the wounded nucleon 
model, based on the idea that the nucleon is composed of two constituents: one 
constituent quark forming a colour triplet and one constituent diquark forming a 
colour antitriplet. We assume, furthermore, (i) that particle production from these 
constituents is independent of the number of interactions they underwent and (ii) 
both constituents produce the same number of particles. It turns out that these two 
assumptions are sufficient to describe correctly the PHOBOS data [1-3] on particle 
production in the central (y = 0) rapidity region. 

In the next section we formulate the model in more detail. In Section 3 the 
determination of the model parameters from proton-proton elastic scattering is de- 
scribed. Particle production in Au-Au collisions is discussed and compared with the 
PHOBOS data in Section 4. Our comments and conclusions are listed in the last 
section. 

2 The model 

As indicated in the introduction, we assume that in the soft collisions nucleon can 
be approximated by a system composed of the two constituents - a quark and a 
diquark - acting independently. 

For the process of particle production we assume that each constituent (quark or 
diquark) which underwent at least one inelastic collision emits a certain amount of 
secondary partons. This number is independent of the number of collisions this con- 
stituent underwent afterwards. As this is obviously a straightforward generalization 
of the wounded nucleon concept we call them wounded constituents. Furthermore, 
we assume that the wounded constituent quark and diquark emit secondary partons 
in approximately the same manner. 

The arguments in favour of such an approach are collected in the last section. 
Here we would like to emphasize, however, that they may apply only to soft colli- 
sions, i.e. those in which the transverse mass of the emitted partons does not exceed 
~ 200 MeV. The analogous limits for the observed hadrons are, naturally, somewhat 
higher. 

It should be also pointed out that the proposed description can only be justified 
far from the fragmentation regions of the projectile and of the target. Therefore in 
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this paper we restrict the comparison with the data to the central rapidity region. 

The consequence of these assumptions, fundamental for our conclusions, is that 
the differential multiplicity of partons emitted in the nucleus-nucleus inelastic colli- 
sions can be represented as: 

dN AB (q+dlj? ( \ , (q+d) 



dy 



= wT a, F + (y) + wr a> F4y), (i) 



where w^ +d ^ and Wg +d ^ are the numbers of wounded constituents (quarks and di- 
quarks) in nucleus A and B, whereas F + , F_ are the differential multiplicities of 
partons emitted by one wounded constituent in A and B, respectively 1 . 

When (1) is applied to nucleon-nucleon collisions we obtain: 

^=»rM+f-wi. (2) 

These two equations summarize the relation between particle production in 
nucleon-nucleon and nucleus-nucleus collisions implied by the model. 

For the symmetric collisions (A = B) we obtain the particularly simple and 
elegant result: 

dN AA /dy _ w { J +d) 

Kaa ~ dN NN /dy~ (3) 

Since the R.H.S. of this equation is independent of the phase-space region where 
Raa is measured 2 , but depends on the centrality of the collision, this is indeed a 
very strong consequence of the model. 

At the vanishing cm. rapidity, (1) and (2) imply a simple relation even for 
asymmetric collisions: 

r> , m _ dN AB (y = Q)/dy + ^ +d) 

tiAB{y = 0) = — — - — — — = — - t— -= — . (4) 



dN NN (y = 0)/dy 2w% +d) 



To make a full use of (3) and (4), it is necessary to evaluate w^ +d ^ and w A q+d \w < " q ; 



B 



as function of the impact parameter of the collision. The corresponding formulae are 
obtained by a straightforward counting of probabilities [4, 8]. For the average number 
of wounded quarks and diquarks in each of the colliding nuclei {w^ +d ^ — W q + Wo) , 
at a fixed impact parameter b, one obtains: 

W ^ h ) = -Aj^ f T ( 6 " s ) f 1 " I 1 " ^ q T aq (s)] A [1 - a ad T ad (s)] A ) d 2 s, (5) 



1 This formula of course is valid for any distribution, not necessarily in rapidity. 
2 Provided it is far enough from the fragmentation regions. 
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where a denotes q or d, <JAA{b) = d 2 <jAA/d 2 b is the inelastic differential AA cross- 
section 3 , T{b) is the nuclear thickness function (normalized to unity). T ab {b) is given 
by: 

T ab (b) = — [ a ab (s)T(b-s)d 2 s, (6) 

&ab J 

where a ab (s) = d 2 a ab /d 2 s are the differential cross-sections of the constituents (in 
impact parameter plane) and a ab = f a ab (s)d 2 s are the corresponding total inelastic 
cross-sections (ab denotes qq,qd or dd). 

Note that the formulae (5) and (6) take into account the impact parameter de- 
pendence of the constituent cross-sections. If this dependence is neglected [a ab (s) = 
a ab 5 2 {s)lT ab {b)=T{b) A 

For the number of wounded constituents in nucleon-nucleon collisions we have, 
similarly, w^ +d ^ = w q + w d with: 

Wq )d = [ h q4 (b)d 2 b, (7) 

&NN J 

where ctnn is the total inelastic nucleon-nucleon cross-section and h a is given by: 

h a {b)= J d 2 Sq d 2 s' q d 2 s d d 2 S ' d D(s q ,s d )D(s' q ,s' d ) 

{1 - [1 - a ad (b + s' d - s a )\ [1 - a aq (b + s' q - s a )]} , (8) 

with D(s q , s d ) being the effective thickness of the nucleon. 

It remains to determine the cross-sections of the constituents and their distri- 
bution inside the nucleon. This demands a detailed analysis of the nucleon-nucleon 
collisions, as discussed in the next section. 



3 Nucleon-nucleon collisions 

The distribution of the constituents inside the nucleon and their cross-sections are 
not known. We propose to determine them from the analysis of the data on proton- 
proton elastic scattering. 

Consider first the inelastic nucleon-nucleon collisions. Following the idea that 
the interaction can be described as the interaction of two independent constituents 

3 For heavy nuclei <JAA(b) = 1, except at very large impact parameters which are of no interest 
even for most peripheral events measured at RHIC. In case of AuAu collisions, using the optical 
approximation, we have verified that <JAuAu(b) = 1 for b < 14 fm, corresponding to W > 5. 

4 We have verified that this is a poor approximation for peripheral collisions. 
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in each of the nucleons, we have [8, 9]: 

1 - a(s q , s d ; s' q , s' d ; b) = [1 - a qq (b + s' q - s q )\ [1 - a qd (b + s' d - s q )\ 

[1 - a dq (b + s' q - s d )][l - a dd (b + s' d - s d )}, (9) 

and 

a(b) = J d 2 s q d 2 s' q d 2 s d d 2 s' d D(s q , s d )D(s' q , s' d )a(s q , s d ; s' q , s' d ; b), (10) 

where s q (s' q ), s d (s' d ) are transverse positions of the quarks and diquarks in the two 
colliding nucleons. 

From the unitarity condition we deduce: 

t el (b) = l- y/l-a(b). (11) 

This allows one to evaluate the elastic and total cross-sections. By comparing them 
with data one can obtain information on the parameters of the model. 

Since the nuclear cross-sections are not sensitive to the exact shape of the impact 
parameter dependence of the constituent cross-sections (as the nuclear radius is 
much larger than that of the nucleon), we parametrized a a b(s) using simple Gaussian 
forms: 

a ab (s) = A ab e~ s2 /^K (12) 

The radii R a b were constrained by the condition R 2 ab = Rl + R 2 where R a denotes 
the quark or diquark's radius (a natural constraint for the Gaussians). 

From (12) we deduce the total inelastic cross sections: a ab = i\A ah R 2 ah and we 
also demand that the ratios of cross-sections satisfy the natural condition: 5 

cr qq : a qd : a dd = 1 : 2 : A, (13) 

expressing the idea that there are twice as many partons in the constituent diquark 
than those in the constituent quark. This allows to express A qd and A dd in terms of 
A qq . 

For the distribution of the constituents we again take a Gaussian: 

D(s q , s d ) = l -^e-^' R2 5 2 {s d + Xs q ). (14) 

The parameter A has the physical meaning of the ratio of the quark and diquark 
masses and satisfies, 1/2 < A = m q /m d < 1 (the delta function guarantees that the 
center-of-mass of the system moves along the straight line) . 

5 We have verified that the main results of this paper are not sensitive to the exact values of 
these ratios. A detailed analysis of elastic pp scattering (and, particularly of this assumption) will 
be given elsewhere [10]. 
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Thus, finally, the model contains 5 free parameters. 

Using this formulation and the formula (11) we have evaluated the elastic and 
total proton-proton cross-sections and adjusted the parameters by demanding that 
(i) total inelastic cross section (ii) slope of the elastic cross section (hi) position of 
the first diffractive minimum in elastic cross section and (iv) height of the second 
maximum in elastic scattering are in agreement with data. 




Figure 1: Comparison of the experimental differential cross section with the results 
of the quark-diquark model. Data at -Js = 30.5 GeV are taken from [11]. 

Data at all ISR energies [11, 12] were analyzed. It turns out that the model works 
very well indeed which is by no means a trivial conclusion 6 . One example of such 
calculation is shown in Fig. 1, where the differential cross section da/dt at the ISR 
energy 30.5 GeV, evaluated from the model, is compared with experimental data 
[11]. One sees a rather impressive agreement. A detailed discussion goes beyond 
the scope of this paper and will be given elsewhere [10]. Here we are interested only 
in the resulting values for <J qq /<7NN and w^ +d ^ which are necessary for evaluation of 
the R.H.S. of (3). 

From the point of view of the present investigation, the most important observa- 
tion is that both the ratio a qq /cr NN and the average number of wounded constituents 
in nucleon-nucleon collisions w^ +d ^ seem almost entirely independent of the details 
of the model (provided that, as explained above, the parameters are adjusted to 

6 For instance, an analogous calculation performed in the model with the assumption that the 
proton consists of three uncorrelated constituent quarks led to negative conclusion [13]. 
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describe correctly the proton-proton elastic data). The obtained values are: 

cr qq /cr NN = 1.147 4- 1.148; w% +d) = w q + w d = 1.182 4 1.186. (15) 

These values 7 , supplemented by the relation (13), are used for evaluation of the 
R.H.S. of (3). 

4 Au+Au collisions 

Having determined the parameters of the model from the elastic pp data, we could 
evaluate its predictions for the particle production in Au-Au collisions which is the 
main goal of this investigation. 




100 200 300 400 

W 



Figure 2: Predictions of the wounded quark-diquark model (for W > 5) compared 
with those from the wounded nucleon model. 

Since the PHOBOS data are presented versus the number of the wounded nu- 
cleons in both colliding nuclei (2w^ = W), we have to calculate also as a 
function of the impact parameter of the collision. This can be obtained from the 
well-known formula [4]: 

w A\b) = -4tt h{b - s) {1 - [1 - <J NN T NN {s)\ A ) d 2 s, (16) 
&aa{o) J 

with (Jnn(s) in a Gaussian form with o-atat(O) = 0.92 taken from our estimates (it 
agrees very well with the data [12]). 

For the nuclear density we have been using the standard Woods-Saxon formula 
with the nuclear radius Ra u — 6.37 fm, and d = 0.54 fm. 

7 Note that <J qq /<JNN > 1/9, indicating presence of shadowing. 
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In Fig. 2 the predicted ratio RauAu/w'au (which shows explicitly the deviation of 



our model from the traditional wounded nucleon model) is plotted versus 2w 



(TV) 
Au 



w 



for W > 5. One sees that the model explains naturally the increase of the production 
multiplicity from one wounded nucleon with increasing centrality of the collision. 

The comparison of the model with the PHOBOS data on particle production per 
one wounded nucleon is shown in Fig. 3. 
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Figure 3: The predictions of the wounded quark-diquark model (for W > 5) com- 
pared with the data from PHOBOS coll. [1-3]. The shaded areas reflect the inac- 
curacy in the pp data. 

The data on particle production in pp collisions, necessary to obtain the model 
predictions, were taken from UA5 collaboration, as quoted in [2,3] 8 . They are also 
shown in the Fig. 3 as points at W = 2. This introduces some uncertainty, as 
indicated by the shadowed areas 9 . 

The inelastic proton-proton cross sections needed for this calculations were taken 
as unn = 32 mb, 36 mb, 41 mb and 42 mb at energies y/s = 19.6, 62.4, 130 and 
200 GeV, respectively. Our model then gives the following values of the ratio of the 
integrated inelastic quark-quark to proton-proton cross sections: cr qq /a NN = 0.147, 
0.148, 0.148, 0.149. Finally, the number of wounded quarks and diquarks in a single 
proton-proton collision w q + Wd= 1.183, 1.185, 1.186, 1.187. 

One sees that, within the experimental accuracies a very good agreement both 
in shape and in absolute value is obtained. 

8 For V« = 62.4 and 200 GeV they are taken directly from UA5 data, for </s = 19.6 and 130 
GeV they are interpolated. 

9 Unfortunately, the pp data from PHOBOS coll. are still not available. 
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5 Comments and discussion 



Several comments are in order. 

(i) The concept of wounded nucleons and of wounded constituents is based on 
two ideas: (a) during the interaction, the "wounded" object acts as one unit and (b) 
particle emission process takes much longer time than it is needed for the projectile 
to pass the internuclear distance. These assumptions can be qualitatively justified 
only for soft collisions where the momentum transfer and transverse mass of the 
emitted partons are small enough. 

The first condition requires that the momentum transferred to the projectile is 
smaller than its inverse size. For the size of order of 1 fm, this limits the momentum 
transfer to about 200 MeV. 

The second condition demands that the emission time: 

t-~ (17) 

where 7 is the Lorentz factor of the emitted particle (parton) in the rest frame of 
the target nucleus and y is its rapidity in the same frame, should be significantly 
greater than the intranuclear distance. For y > 2 this limits the transverse mass of 
the emitted partons below 200 MeV. Of course for the observed final hadrons this 
limit may be significantly higher. 

These estimates are, surely, rather crude. A more detailed verification of the 
model for particles with various masses and transverse momenta will be of great 
interest, as it may help to understand better the very concept of a wounded con- 
stituent. 

(ii) As seen from these arguments, the model is not expected to apply in the 
fragmentation region of the projectile and target where, moreover, important energy- 
momentum conservation effects, as well as secondary interactions inside the nucleus 
must be present. Therefore we focus our attention on the central rapidity region 
which, at RHIC energies, is well separated from the fragmentation regions. 

(iii) The model assumes that a wounded quark produces the same number of 
secondary partons as a wounded diquark. This is not unreasonable since the colour 
content of both constituents is the same (3 and 3) and that, probably, the colour 
charge of the projectile is the main factor determining the emission intensity. We 
admit, however, that since the actual dynamics of the soft production process is not 
yet understood, this argument can be questioned. The good agreement of the model 
with data, as presented in this paper, may thus serve as an (indirect) confirmation 
of the important role of colour dynamics in the process of particle production. 

(iv) The predictions of the model described in this paper refer to the emission of 
"primary" partons and do not take into account further development of the system 
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during its expansion and final formation of hadrons. The observed agreement with 
data indicates that the space-time development of the system, despite presence of 
the well-known collective effects, does not introduce drastic changes in its global 
characteristics. This may be considered as an argument in favour of the laminar 
hydrodynamic expansion suggested already by MC simulations of this process which 
seems to indicate very small viscosity of the created medium [14]. 10 

(v) We have verified that the main properties of the quark-diquark picture of the 
nucleon obtained from analysis of the pp elastic data are not sensitive to the details 
of the calculations (several forms of the distribution of constituents and of their 
cross-sections were considered). The typical values of the parameters are R q ~ 0.3 
fm, R d ~ 0.75 fm, R w 0.3 fm, A dd w 0.55, A qd ps 0.5, A qq « l. u Thus in our model 
the diquark appears to be rather large, comparable to the size of the proton. It is 
remarkable that this feature agrees well with other estimates [15], based on rather 
different arguments. 

(vi) In this paper we have only discussed the symmetric Au-Au collisions. It 
would be also interesting to check the model for asymmetric collisions. We have 
verified that in case of d-Au collisions at y = the formula (4) gives the result 
which does not differ very much (less than 10 %) from that of the wounded nucleon 
model. Thus the observed good agreement of the wounded nucleon model [16] with 
the PHOBOS data on d-Au collisions [17] is not destroyed in our approach. 

In conclusion, we have formulated a model in which the soft collisions of the 
nucleon are described in terms of interactions of its two constituents: a quark and 
a diquark. The model can be adjusted to describe very precisely the elastic proton- 
proton scattering data. Supplemented with the idea of wounded constituents, the 
model accounts rather well for the centrality dependence of particle production in 
the central rapidity region at RHIC energies. 
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